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The problem of correlation analysis in MST radar is to determine the 
scattered power, Doppler frequency and correlation time for a noisy signal. It 
is assumed that coherent detection has been employed, with two ? accurately 
balanced quadrature receiving channels. It is further assumed that coherent 
integration has been performed with a window length significantly less than the 
correlation time of the signal. 

The analysis problem may be looked at from the point of view either of 
Fourier analysis or of correlation analysis, and it must be emphasized that the 
two approaches, if used properly, give identical results. Why. then, use 
correlation analysis at all? The reason can be seen from the spectrum and 
correlation function shown in Figure 1. In each case 1 min of data is 
represented, for example (with 1/8 sec coherent integration time) 480 pairs of 
data points (real and imaginary). Ordinary discrete Fourier analysis requires 
about 2 x 10^ floating-point multiplications, all involving transcendentals. 

In fact, however, only the area P under the echo spectrum, its position f^ 
and its width f 2 are required; these quantities have to be calculated from the 
spectrum by separate algorithms. 

It can easily be shown that the quantities P, fj and f 2 can be determined 
from the first few values of the complex autocovariance functions, shown on the 
lower part of Figure 1. This function can be calculated out to a number of 
lags approaching the length of the sample, but almost no additional information 
is contained in the part of the function beyond the first few lags. In 
examining such a function, it is necessary to make an assumption; namely, that 
the curve for lags other than zero can be extrapolated back to zero to give the 
signal power P and the noise power N as shown; this is possible because the 
noise power is uncorrelated from one coherently integrated sample to the next, 
while the coherent integration time has, as indicated above, been chosen so as 
to make the correlation between one sample and the next very good. 

The correlation time (or time to correlation equals .5) is estimated as 
shown after the noise power has been removed. The spectral width f^ is the 
reciprocal of the correlation time. The Doppler frequency f^ is found from the 
slope of the imaginary part of the correlation function at the origin, care 
being taken to eliminate the unwanted variance N. 

In principle* only 3 points on the complex autocorrelation function need to 
be calculated, which would require about 6000 multiplications, many fewer than 
the number required for the spectral approach. In fact, a larger number of lags 
(up to 12) is often calculated with the idea of improving the analysis. As 
described by COUNTRYMAN and BOWHILL (1979), the values of the argument of the 
complex covariance for the various lags can be weighted to give a more accurate 
value. A maximum likelihood analysis of the estimation problem is given in 
Appendix 1. This method of analysis has not yet been applied to experimental 
data. 


The calculations of autocorrelation functions, , however , need not 
necessarily use multiplication at all. BOWHILL (1955) described a method of 
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finding autocorrelations using the mean difference of separate samples rather 
than their mean product. Appendix 2 gives a description of the algorithms used. 
This technique has been successfully applied in a microcomputer operating system 
for the real-time processing of MST data (see paper 8 .3-D). 

There are other ways in which the correlation process can be speeded up. 
HAGEN and FARLEY (1973) describe 11 methods. Table 1 illustrates several of 
these algor ithms, together with the efficiency in terms of use of the input 
data. It should be emphasized that the technique of Appendix 2 has an 
efficiency in excess of 95%. The question arises as to whether a hardware 
correlator is worthwhile for MST radar* It is my opinion that the use of 
coherent integration, which reduces the number of input data by about 2 orders 
of magnitude, makes a special hardware device unnecessary , particularly with the 
use of algorithms such as those described in this paper. 




Figure 1. Power spectrum and autocovariance function 
corresponding to a coherently scattered signal plus 
noise. 


Table 1 

Type of Correlation 

Output 

Efficiency (%) 

Multibit 

■® 2 .p 

100 

One-bit x multi bit 

(2/ir)' L ^ 2 op 

64 

One-bit 

(2/tt) sirT^p 

41 
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APPENDIX 
Suppose that 


1. Maximum-Likelihood Estimation of C-S Parameters. 

the real part of the complex autocorrelation function is given 


by 





r(x) «= r Q exp(- 

■ax 2 ) 

costox 

where 

e 0 (T) = l-r 0 

(t 

* 0) 


= 0 

(t 

> 0) 

and 

t “ 0, 1, 2 .. 

... n. 



Now let the observed autocorrelation function have real part R(x) . For the 
imaginary part, i(x) « r 0 exp (-ax 2 ) sinwx + e 0 (x) if the r and i channels have 
equal sensitivity. 


Now let the noise n on the r and i channels have a Gaussian distribution 
exp (-bn 2 ) then the likelihood of a given set of l(x) and R(x) is 


and 


L - exp [ -b t (R(t ) - r(r)} 2 - b E {I(t) -i(x)}] 
o o 

-knL»Z [R(t ) -r(x)l 2 + I [X(t) -i(x)] 2 
b o o 


and maximum likelihood amounts to a least squares fit of RCO and I(x) by r( T ) 
and i(x). 

Substituting for r(x) and i(x), and neglecting x - 0, 

- tr-ln L = £ [R(x) -r exp(-ax 2 ) coswt ] 2 + l l(x) -r exp(-ax 2 ) sinwx] 2 
d ^ o j o 

* £ [R 2 (x ) + I 2 (x) + r^ 2 exp(-2ax 2 ) - 2r Q 2 exp(-ax 2 ) (R(x) cosoix + l(x) simix)] 
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Differentiating with respect to r Q , a and to, and equating to 0, we get the 
following equations which must be satisfied simultaneously: 

n on _ 

r £ exp(~2ax z ) - £ exp(-ax z ) [R(x) coaox + I(t) sintox] = 0 


or 


r n £ x 2 exp(-2ax 2 ) - I T 2 exp(-ax 2 ) [R(x) coswx + X(x) sintox] *= 0 
1 1 


£ exp (-ax 2 ) [R(x ) sintox - l(x) costox] * 0 


LIhl for UT « i. 

£ R(x)x exp(-ax z ) 


APPENDIX 2. Rapid Pseudocorrelation Technique. 


Consider two voltages R(t) and I(t), nominally the real and imaginary parts 
of the phasor of a radar return signal. The problem is to determine the center 
frequency and power, and correlation time, of an embedded signal of frequency 

REPRESENTATION OF R(t) AND l(t) 

The signal of frequency to can be represented by a phasor S(x) at that 
frequency, giving S^t) coscot and S-r(t) cos(tot+4>) in the voltages R(t) and 
I(t), respectively. S^(t) and Sj(t) are random variables, such that 
S^CtJ/SxCt) « constant. The phase shift <j> is nominally ir/2. A noise voltage 
n(t) will also appear in R(t) and I(t), which is supposed to be completely 
uncorrelated from one pulse to the next. 

We therefore have 

R(t) * %(t) co scot + n R (t) 

I(t) * Sx-(t) cos(tot+(|>) + n-j-(t) 

MEAN SQUARE DIFFERENCE DEFINITIONS 

Let RO 2 
I0 2 
RR 2 (t ) 

II 2 (t ) 

RX 2 (t ) 

IR 2 (t ) 


- <R 2 ( t ) > 

= <X 2 (t)> 

“ <[R(t) - R(t + t)] 2 > 

- <[l(t) - I(t + t)] 2 > 
m <[R(t) - l(t + T )] 2 > 

- <[I(t) - R(t + x)]2> 
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Now all these quantities may be related to the mean absolute values of the 
squared quantities by the relation: 

<X 2 (t) > - k[<|x(t) |>] 2 

where k is a constant for a given waveform. 

EVALUATION OF DIFFERENCES 

From the expression for R(t ) , 


<R 2 (t) > = -| <S R 2 (t) > + <n R 2 (t) > 

since s R (t) and n R (t) are independent random variables. Further, we define 


s R 2 = <S R 2 (t», N r 2 = <u R 2 (t) >, 

RO 2 = § S R 2 + N r 2 
I0 2 = \ S R 2 + N R 2 

2 

RR 2 (t) = <[S R (t) coswt + n R (t) - S R (t + t) cos(wt + wt) - n R (t + r)] > 

= ~ S R 2 + N r 2 + ~ S R 2 + % 2 - 2<S R (t ) S R (t + t) coswt cos(wt + wt)> 
= S R 2 + 2N r 2 - <S R (t) S r (t)> co S oot 


and defining 


p ( t ) = <S R (t) S R (t + T )>/ S R 
RR 2 (t) = S r 2 [1 - p (t ) costur] + 2N r 2 
II 2 (t) » S^U - p(t) cosojt] + 2N-J- 2 

Similarly , 

o 2 

RI (t) = <[S R (t) co 8io t + n R (t ) - S^(t + x) cos(u>t + wt + <j>) - Hj(t+T')] > 

1 ? 7 12 2 

= 2 S R + N R + 2 S I + N I 

+ 2 <S R (t) Sj(t + t) co scat cos(wt + wt + 0 )> 

+ 2 <n R (t) n^t + t)> 

and the latter term is always zero. So 

7 1 2 1 9 O 9 

RI (t) = j S R + 2 s i + n r + N i 

+ 2 <S R (t) S R (t + t) co 8 (tot + <(>)> 

= -k S 2 + -I S 2 + N 2 + N 2 + S n S T P(x) cos(+o>t + <)>) 

2 R 2 1 R 1 R 1 


Similarly, 


IR 


2 (t) = \ S R 2 + \ S x 2 + N r 2 + Nj 2 + SpSj p(t) cos(-u)T + ()>) 
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POWER AND FREQUENCY CALCULATIONS 

From the above relations, 

2R0 2 - RR 2 = S R 2 p(t) co but 
2 2 2 

2X0 - II - Sj p (x) co suit 

RI 2 - RO 2 - IO 2 * S R Sx P(t) cos(wt + 4>) 

IR 2 - RO 2 - IO 2 = Sj^Sj p(t) co s ( -ibt + f) 

and 

I ( 2R0 2 - RR 2 )(2I0 2 - II 2 ) ] 1/2 - Sj^SjPCt) <®but 

IR 2 - RI 2 = 2 SgSjP (t ) sinoT • sinsj) 

RI 2 + IR 2 - 2(R0 2 + IO 2 ) - 2 S r S i P(t) coswt • cos4> 

Normally, $ is set to it/2, so the first two equations become 

[ (2R0 2 - RR 2 ) (210 2 - II 2 )] 1 / 2 - S D S T p{x) COS60T 

K 1 

IR 2 - RI 2 = 2 S p S t p(t) sinuT 

and and wt may be found trigonometrically. S R S j may be adopted as the 

measure of scattered power. The third equation may be used as a check upon the 
accuracy with which the phase-quadrature channels have been adjusted. 



